Dissipative and stochastic geometric phase of a qubit within a canonical Langevin 

framework 
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Dissipative and stochastic effects in the geometric phase of a qubit are taken into account using 
a geometrical description of the corresponding open-system dynamics within a canonical Langevin 
framework based on a Caldeira-Leggett like Hamiltonian. By extending the Hopf fibration S A — > 
S 2 to include such effects, the exact geometric phase for a dissipative qubit is obtained, whereas 
numerical calculations are used to include finite temperature effects on it. 



Introduction. The concept of geometric phase (GP) 
in quantum systems was proposed by Berry [l[ when he 
studied the dynamics of an isolated quantum system that 
undergoes an adiabatic cyclic evolution. This cyclic evo- 
lution is due to the variation of parameters of the Hamil- 
tonian and is accompanied by a change of the wavefunc- 
tion by an additional phase factor which depends only 
on the geometric structure of the space of parameters. 
The underlying mathematical structure behind GPs was 
pointed out almost simultaneously by Simon Q. Soon 
after, the generalization to non-adiabatic cyclic evolu- 
tion was carried out by Aharonov and Anandan Q and 
non-cyclic evolution and sequential projection measure- 
ments by Samuel and Bhandary Q . Although GPs have 
been observed in the laboratory [H-0]' realistic quan- 
tum systems are always subject to decoherence due to 
their surroundings. Therefore, the extension of the GP 
to the case of open quantum systems becomes funda- 
mental. The first formal extension of the GP was car- 
ried out by introducing the concept of parallel transport 
along density operators Q. In a more physical context, 
the concept of GP was generalized for non-degenerated 
mixed states and for degenerate mixed-states under 
unitary evolution [l(| • Using a kinematic approach fill ] . 
GPs for mixed states in non-unitary evolution were ad- 
dressed [lH. Within a spin-boson model, GPs in open 
quantum systems have also been calculat ed [l3i together 
with a study of their geometric nature [lj]- A differ- 
ent approach was introduced in Ref . [IH , where the GP 
was described by a distribution. In classical physics, the 
counterpart of the Aharonov-Anandan (or Berry) phase 
was early discovered by Hannay [l6| . Regarding classical 
dissipative systems, GP shifts have been defined in dis- 
sipative limit cycle evolution [l7l - IT9| . showing that they 
can be identified with the classical Hannay angle in an 
extended phase space [20|. Moreover, GPs can be con- 
structed for a purely classical adiabatically slowly driven 
stochastic dynamics (2lT[23| . 

In this Letter we tackle the problem of including dis- 
sipative and stochastic effects in the GP of a qubit. Our 
study will be based on a geometrical description for a 
non-isolated qubit within a canonical Langevin frame- 
work (see [24| and references therein) using a Caldeira- 
Leggett like Hamiltonian (25| . 

Mathematical preliminaries. It is well known that if 



|\&) represents a normalized n-level system, then 1$) € 
g2n-i Thus, the geometry of odd-dimensional spheres 
is related to the quantum mechanics of finite-dimensional 
Hilbert spaces. In fact, the celebrated Hopf fibration [26[ 
relates the quantum and classical description of qubits 
by means of the map ir : S 3 — > S 2 . This map can 
be understood as a composition, tt = 3 o 0, where 
Q : S* 3 C C 2 CP 1 sends an element of C 2 to its equiv- 
alence class and 5 : CP 1 (= CUoo) — > S 2 is given by 
the stereographic projection. It can be shown that the 
Hopf map can be written in terms of the Pauli matrices 

as tt (i*) G s 3 ) = m^m, (*|*«|*>, (*|^|*» € S 2 , 

where (^Ict^) 2 + (^\a y \^) 2 + (*|<r 2 |*) 2 = 1. Thus, 
from the Hopf map it can be shown that quantum and 
classical mechanics may be embedded in the same formu- 
lation. Specifically, for a qubit, the Strocchi map [27[ is 
exactly the Hopf map previously described. However, the 
S 7 — > S 4 Hopf map, which is an entanglement-sensitive 
fibration, does not have a classical analog (2§J. As there 
is a map S 12 ™ -1 — > CP™ -1 and the complex projective 
space has a natural symplectic structure (CP" -1 is a 
Kahler manifold), n-state systems have a well-defined 
classical correspondence, which is given by the Strocchi 
map [27j |. Thus, one can derive a classical Hamiltonian 
function for a n-level system by defining, for example, 
appropriate n action-angle coordinates in CP™ -1 [291 ] . 
For example, as shown in (30j . the symplectic structure 
of S 2 is responsible for the Aharonov-Anandan GP. 

Using the pair of action-angle coordinates (I, <E>) on 
S 2 , the Hopf map can be expressed as ^(1^)) = 

ai\l) + 



s/T- 



I 2 sin$, /) 



where 

|2 



and $ = 



- I 2 cos$, 

a 2 \2) (aj = |%|e^ G C), l = \at\~ - \a 2 \ 
4>i — 4>2- Thus, the Hamiltonian operator H = J^iVi&ii 
where &i are the Pauli matrices and r\i € R, can be 
Hopf-mapped to a Hamiltonian function, Ho, given by 
Ho = -2\/l - I 2 (jn cos $ + 772 sin $) + 2p 3 I (this is the 
Meyer-Miller-Stock- Thoss Hamiltonian [33, [32[, widely 
used in molecular physics). 

Dissipative-stochastic Hopf fibration of S 3 . Our study 
will be based on a Caldeira-Leggett [231 like Hamilto- 



nian for a qubit in the Langevin framework (see }24| and 
references therein), which can be expressed as 



c z Xi+^2<5> 2 c 2 



(1) 



2 



where the oscillator mass has been taken to be one and 
Cj are the system-bath coupling constants. This model 
takes into account a renormalization term due to the in- 
teraction with the environment. 

Let us start with pure and normalized quantum states. 
Using the (J, 4>) action-angle coordinates, it is easy to see 
that the unit-radius sphere S 2 is the set of points satis- 
fying (-■ v/l -P cos$) 2 + (-Vl - I 2 sin$) 2 + I 2 = 1. 
The radius of this sphere remains constant along the 
time due to energy conservation. That is, the cor- 
responding Hamiltonian function, given by Hq = 
— I 2 (cos $ + sin $) + I remains constant in time 
(we have taken rji = 1 for simplicity). When Ohmic 



dissipation is assumed, and after removing the bath vari- 
ables, the corresponding equations of motion issued from 
Eq. (P) are 



/ = - s/l - I 2 (sin $ + cos $) - 2 7 $(f ) + £(t) 
I 



$ = 



(cos<I> + sin$) + 1, 



(2) 



which correspond to the following effective Hamiltonian 
function H(i) = - \/\ - P (cos $ + sin $) + I + 27$$ - 
£(£)$. This stochastic dynamics can be interpreted in 
terms of a stochastic Bloch sphere /f> 2 (7, £) defined by 



yl - P cos$) + ( VI - J 2 sin $ 



i 2 = i-i Jt ^ 2 )+am 



(3) 



where 7 is the friction constant, is a stochastic Gaus- 
sian process. The time-dependent radius is given by 

R t(liQ = 1 - lf t (® 2 ) + Thus > dissipative and 

stochastic effects make the Bloch sphere breathe, Eq. (|3|) , 
by changing its radius in time. This radius is bounded 
for £ = 0, Rti'Y, 0) < 1. The equality is reached at t = 
and at asymptotic times, where thermal equilibrium is 
reached (in this case, a point in the unit radius sphere 
represents a pure state). Otherwise, mixed states are 
represented at each instant of time as points in different 
spheres of variable radius, as shown in Fig. ([1]). When 




Figure 1. Representation of the dissipative-induced time de- 
pendence of the Bloch sphere for £ = 0. Time is running in 
the radial direction. 

I 



£ 7^ 0, the square radius becomes a stochastic variable 
which could reach i? t ( 7 ,£) > 1 for a particular phase- 
space trajectory. 

The breathing of the Bloch sphere can be geometrized 
by extending its round metric by adding both dissi- 
pation and noise. If these terms are included, the 
metric of S 2 can be written as ds 2 [-S 2 ( 7 , £)] = 

R 2 t {j,0ds 2 [S 2 ], where ds 2 [S 2 ] = ^ + (1 - P)d<S> 2 
is the round-metric for S 2 in action-angle coordinates. 

In order to calculate the stochastic GP, let us start 
with pure states. If we choose an orthonormal mov- 
ing frame field er 1 



dl 



2 = VI - Pd$, then 

da 2 = — Tr-fjdl A # = -la 1 A d$ = Id<S> A a 1 . Us- 

v l— 1 

ing Cartan's first structure equation, we obtain the only 
non-vanishing connection one-form, oj = /e?$. Thus, the 
dynamic phase is obtained as 4>d = § u> = § Id<f> = irl = 
7rcos0. If S 2 is considered as an embedded submanifold 
of S" 3 , taking into account that the spin connection of 
S 3 is £0 = dE + Id&, where S is the extra Euler angle 
which parametrizes the third dimension, it can be shown 
[28j that S 2 inherits from S 3 the connection one-form 
Q = (I — l)d$>. Thus, the GP is expressed as 4> g = §Cj. 
Now, if dissipation and noise are taken into account, the 
orthonormal moving frame field is given by -fft( 7 ,£)(7 l 
with i = 1,2. In this case, the only non- vanishing con- 
nection one-form is a; (7, £) = R t (j,£,)Id& = i? t (7,£)w. 
This allows us to define the stochastic dynamic phase ac- 
quired after a cycle of period T as 4>d{li £) = § w (7i £) = 

§R t { 1 ,ti)I{T)d$> = cos0(r)/ v /l- 7 |($ 2 )+^d$. 
Therefore, the stochastic connection one-form, which can 
be defined as w( 7 ,£) = [Ji?t( 7 , £) — 1] d^, leads to the 
corresponding GP 



[I(T)Rt{7,0-i\ 



HT)\1 



- 1 



(4) 
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Notice that the dynamics driven by this kind of Caldeira- 
Leggett like coupling could be interpreted in terms 
of a mapping between conformal spheres with con- 
formal unitary group fibers at each instant. More- 
over, the GP defined by Eq. ((4]) becomes </> 3 (7,£) = 

I J 1 - 7^-(* 2 ) + - 1 ] d<& under the transforma- 



is a 



tion I*) = e m |*) (or $ = $ + a). Thus, ff (7,£) 
gauge invariant quantity. 

Dissipative qubit (zero temperature). Let us consider 
a simple qubit which can be represented by the Hamil- 
tonian operator H = ea z . The corresponding dissipative 
dimensionless Hamiltonian function (t — > 2et) , which can 
be written as H(t) = I + i^.$>&, leads to the pair of 
coupled equations I = — ^$ and $ = 1 with solutions 
I(t) = I - 3-t and <&(t) = $o + *, where Iq and <J> are 
the initial conditions of the action-angle variables. After 
a cycle of period T — 2ir (the scaled Rabi frequency is 
luq = 1), the dissipative GP can be computed as 



J fr) = f[Rtb)I-l]d$ 



I(T) 



1 - — <M> 



= -7T + - 



2e 

TTJ\ 3 / 2 

27; 



7T COS Oq 

7 



(5) 



where I = cos 9 . Notice that the non-dissipative GP is 
recovered when 7 — > 0. 

It is also interesting to remark that Eq. ([5} is only valid 
for 77r/2e < 1, which reflects nothing but the Caldeira- 
Leggett energy rcnormalization in an Ohmic environ- 
ment. Moreover, this energy renormalization can be re- 
interpreted in terms of a Block sphere- whose radius de- 
velopes harmonic oscillations. In order to show this state- 
ment in simple terms, let us introduce the damping fac- 
tor by means of the phenomenological Caldirola-Kanai 

I 



Hamiltonian [33l435| for a harmonic oscillator (which is 
equivalent to the corresponding Caldcira-Lcggctt Hamil- 
tonian for the oscillator at zero temperature). The 



Hamiltonian reads H = (l/2)p e 



9 

+ q e^ 



(the 7r 



factor has been introduced to note that time evolution is 
cyclic) . It is straightforward to derive the corresponding 

equations of motion, leading to y + I — (-57) 2 y = 0, 

where y = p,q. Thus, the renormalized frequency due to 



the damping term is ui 2 



which has physical 



sense only when 77r/2e < 1, as Eq. ([5]) shows. There- 
fore, the radius of the Bloch sphere developes harmonic 
oscillations with a renormalized frequency (another in- 
terpretation is that < 1 requires 77r/2e < 1). 

The weak coupling limit of the GP given by Eq. ([5]) 
can be expressed as 



<t>g{l) = -71" (1 



cos*b)-2(f 



(cos #0 



4) + 0( 7 2 ). 

(6) 

A direct comparison with other authors [10|, [HI, [l5| is not 
pertinent since different models of dissipation and other 
approaches to calculate the GP were used. 

On the other hand, for this dissipative dynamics, in- 
formation on interference experiments can be straight- 
forwardly extracted from the probability density itself. 
The typical interference intensity, J t , evolves in time 



according to J t oc ^(i)^ 



1 + y/l - P(t) cos$(i) = 

1 + a/1- (J - £t) 2 cos(t + $o)j depending critically on 
the ratio 7/2e which is directly related to the GP given 
by Eq. ©. 

Stochastic qubit (non-zero temperature). Noise ef- 
fects can be included in a simple way by assuming a 
Gaussian stochastic process with distribution function 

p(£]0) — \/Jf CX P (~^2~) ' WnCrC P * S * nC mverse 01 

the temperature, which is given in units of 2e. Thus, 
for a qubit, the squared radius is also a stochastic pro- 
cess given by i? 2 (7,£) = 1 + (£ — ^)$. Therefore, the 
stochastic GP acquired after a cycle of period T is given 
by 



(7) 



r 



The previous integral in the noise variable has no analytic solution. Thus, we firstly integrate in $ obtaining 
I 



tt. (8) 



r 



To illustrate the computation of this integral (which docs not have any analytical solution either) , only the GP for 



4 



a non-dissipative qubit at finite temperature will be con- 
sidered. Moreover, as the noise term becomes more im- 



portant for high temperatures, we assume that the mean 
of the Gaussian process is /3 _1 . Therefore, the GP can 
be factorizcd as 



(j) g (/3) = cos6> 



2 

3? 



3/2 



caa8 f(p)-ir, 



(9) 



where /(/?) encodes the thermal information of the GP 
acquired by an isolated qubit. This function is depicted 
in Fig. ([2]) within a large range of temperatures, dis- 
playing a linear behavior at high temperatures and a 
crossover at a critical temperature given by T c « 2e 
(T c pa 1 if adimcnsional temperatures are used). This 
temperature corresponds to the energy difference of the 
two levels of the qubit. Moreover, as f{(3) — > 7r at very 
low temperatures, the GP acquired for the pure state at 
zero temperature is recovered. Finally, notice that this 
type of calculations can be straightforwardly extended to 
any dissipation value. 




1e-07 1e-06 1e-05 0.0001 0.001 0.01 
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Figure 2. Thermal behavior of the GP obtained by numeri- 
cal integration of Eq. ([9]) (double logarithmic scale is used). 
Notice that f(/3) — > n at very low temperatures. 



In summary, a geometrical description of the dissipa- 
tive and stochastic dynamics of a qubit within a Langevin 
formalism (with a Caldeira-Leggett like coupling) has 
been developed. The Hopf hbration S 3 — > S 2 has been 
extended to include both stochastic and dissipative ef- 
fects in terms of a Bloch sphere which develops har- 
monic oscillations (in the dissipative case). This proce- 
dure has allowed us to define a gauge invariant geometric 
phase, which has been computed for both dissipative and 
stochastic cases. 
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